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The one-boson exchange potential model for Ac A*' is constructed and the possibility of A.cN bound 
states is examined. We consider an effective Lagrangian for the charmed baryons, Ac, Ec, and Ec, 
reflecting the heavy quark symmetry, chiral symmetry, and hidden local symmetry. We determine 
the coupling constants using various methods. With the derived nonrelativistic potentials, we study 
the bound state problem of the AcA'^ system and relevant coupled channel effects. It is found that 
molecular bound states of A.cN are plausible, for which the channel couplings of ScA^ and EcA^ are 
essential. 
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I. INTRODUCTION 



Heavy quark hadrons attract great interests since the observation of the exotic Z?sj(2317) and X(3872). One of 
their prominent features is the proximity to one or several thresholds. These mesons, followed by the observation of 
other exotica [iH^i triggered lots of discussions about the bound state problem of heavy quark hadrons. The case of 
' (-H hidden charm mesons is a major concern because they are relatively easy to produce in the experimental facilities. 
O There are several types of two-body hidden charm systems: [cq)(cq), {cqq){cqq), {cq){cqq), {cc)(qq), and {cc){qqq). 

I For the first three types, the components are heavy and may be treated nonrelativistically. Since the interactions 
may be mediated by light meson exchanges, each type had been proposed to have bound states 0-[l3|. For the last 
(— I . two types, the two-body attraction is not so strong and the formation of a molecule may be difficult. However, it is 
I— I ' possible that the in- medium effects may drive the formation of bound states because of the van der Waals type force. 
Such quarkonia are called nuclear-bound charmonium [l8| or hadro-charmonium [loj . 

More exotic bound states are possible in the heavy quark case due to the relatively small kinetic term. The charmed 
^ ^ meson-nucleus systems are interesting in that they are easy to tag . One may find works about DN (or DN) interaction 
• and the charmed mesic nuclei in various methods in Refs. |20l429l |. The calculations favor the existence of such bound 
\^ ' states. 

■rj" . The formation of the Ac and Ec hypernuclei was predicted more than 30 years ago |30|. There were further studies 
' about such hypernuclei [3ll - [33l |. The observation of Ac-nuclei had also been claimed [34|. Recently, this problem 
, was revisited with the quark-meson coupling model [35| and it was found that the charmed hypernuclei are still 
T-H ■ probable. However, it is inconclusive whether the two-body bound states in the AciV channel exist or not (soj . In 
Refs. m, [13, various heavy quark deuteronlike bound states, A^Sc, A^S^, iVScc, SScc, and so on, were discussed 
by scaling the strengths of the nucleon-nucleon interaction models. But the interactions between Ac and N could 
not be approximated by the rescaled A'^A'^ interaction. As AcA'^ is the most fundamental charmed nucleon state, it is 
/\ worthwhile to give it a serious study with the modern heavy quark effective theory and channel couplings to TicN and 

- - - With the forthcoming facilities such as GSTFAIR and J- PARC [38[, it would be possible to draw a conclusion 
whether the charmed hypernuclei exist or not. To understand further this problem theoretically, we revisit the two- 
body bound states in the AcA^ channel at the hadron level in this paper. Since the one-pion-exchange interaction 
is forbidden in the AcA^ channel, we consider scalar and vector meson exchange contributions. We also include the 
coupled channel effects due to the nearby 'ScN and EcA^ thresholds. The coupled channel effects in the heavy quark 
scenario probably have sizable contributions even though the thresholds are not very close to each other, which may be 
realized from the discussions about the resonance effects for the DD production [33 E^- It is interesting to investigate 
how important the coupled channel effects are for the formation of the AcA'' bound states. 

In the early studies of the AcA^ interaction [s^l, the SU(4) symmetry was used to extend the one-boson-exchange 
model for A^A^ and YN interactions. At that time, the heavy quark symmetry had not yet been explored. With 
this symmetry, one does not need to consider the largely broken SU(4) any more. As a result of the heavy quark 
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symmetry, parts of interacting terms are related and the Lagrangian can be written in a compact form. Because 
of the spontaneous symmetry breaking of QCD, it is natural to consider chiral symmetry for pions. To include the 
light vector mesons, we turn to the hidden symmetry approach where the vector mesons are treated as gauge fields 
of a hidden local symmetry. This approach has been successful in explaining the KSRF relation = 2/^gp^^, the 
universality of p coupling, and p dominance of the photon coupling to pions |41[. Here, when deriving the one-boson- 
exchange potentials, we use effective Lagrangians satisfying the heavy quark symmetry for heavy baryons, chiral 
symmetry for pions, and hidden local symmetry for vector mesons. The Lagrangian is analogous to the heavy meson 
version presented in Ref. [4^ . As for the unknown coupling constants, we turn to the quark model, chiral multiplets, 
vector meson dominance, and QCD sum rule. 

Extension of the study to the bottom case is also interesting, but we would not explore it here. Our paper is 
organized as follows. After the introduction, we present the effective Lagrangians in Sec. [Ill Then we determine the 
coupling constants with various methods in Sec. IIIIl In Sec. IIV[ we present the derived potentials. The numerical 
results for the 5- wave spin-singlet and triplet KcN are given in Sec. |V]and Sec. I VII respectively. Finally, we present 
our discussions and conclusions in Sec. IVIII 



II. THE LAGRANGIAN 



Here, we consider the S'-wave KcN states and the channels which couple to them, i.e., the I = \ and = 0+ or 
1+ two-baryon states. We consider the contributions from the YicN and S]*iV channels. It is a 3-channel problem for 
= O"*" and 7-channel problem for = 1+. The labels are listed in Table H] Because of the higher mass, we assume 
that the A(1232) contributions are negligible. 



Channels 1 2 3 4 5 


6 7 


= 0+ KcN{^So) ScNCSo) E*N{''Do) 

= 1+ AcNi^Si) Y;^N(^Si) KN^Si) AcNCDi) J^cNCD^) 





TABLE I: The S-wave AcN states and the channels which couple to them. 



In the heavy quark limit, the ground state heavy baryons Qqq form SU{3) antitriplet with J 



degenerate sextets with = (i, |) + . As in Ref. [43, we use B^, Bq, 
down the compact form of the Lagrangians, we use the notation of the supcrficld S^, which is defined by 



^ and two 

and Bq to denote these multiplets. To write 



S^^B;^+6^{^^+v^h'Be 



(1) 



where w^j is the 4- velocity of the heavy baryon and S is an arbitrary phase factor. One finds that an appropriate choice 
of (5 is —1 (see Appendix [K^ . which is different from the one often used in the literature. Actually, this phase does 
not affect the final results because it appears only in the transition potentials. In the multichannel case, it is easy to 
prove that the binding energy does not change with the following replacement for the Hermitian Hamiltonian 
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(2) 



where the subscripts of H denote the channels and Si (i=2,3,...) are arbitrary phases. Because of the same reason, 
the convention of the relative phase between the sextet and the antitriplet will also not matter. 

One constructs the effective Lagrangian according to the chiral symmetry, heavy quark symmetry, and hidden local 
symmetry (ill, li^, 0, EH . The coupling terms are constrained by the couplings at the quark level: 0+ — 0+ — DJl, 
1+ — 1+ — 971, and 1+ — 0+ — COT. Here 0+ and 1+ are the spin-parity quantum numbers of the light diquark inside the 
baryons and 5Jt represents the pseudoscalar, scalar, or vector meson. The constructed heavy quark baryon Lagrangian 
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t^B — -^Bj + + ^irit , (3) 

I^B, = itr[B3(ii>.i^)i?3]+i/3str[%>^(r^-y^)B3]+^Btr[B3ai?3] (4) 

+^/35tr[5^«a(^" - V'^)S^\ + AstriS^F^'^S,] + ^5tr[5^a5^] (5) 
Ant = .94tr[5^A^B3]+iA/e'"'^''t;^tr[5,^^AKB3] + /i-C-, (6) 

where we use the notations, Eqs. ([7))- ((T0)) . As = A/g — Afg is the mass difference between the sextet and the antitriplet, 

and CT is the scalar singlet meson. Note that Be (5 ) and Bg (1 ) are degenerate because of the heavy quark spin 
symmetry in this effective theory. In the coupled channel calculation of the AcTV — T^^N — YT^N system, we will use 
the empirical values of the masses of Ac, Sc and E* and therefore the Ec — E* splitting is properly included. The 
coupling constants g\ and 174 are the same as those in Ref. (43| . Other definitions are given below. 

y++ J_y+ J_^'+ \ / J_v*+ 1 "*+ 




We use / = 92.3 MeV for the pion decay constant. The constant gy = rapj {spl^T^) = 5.8 is derived with the 
vector meson dominance (VMD) (ill, IH, liji ■ For the nucleon-nucleon interaction part, we use the following SU(2) 
Lagrangian, 



where is the Pauli matrix, representing the isospin. 



III. THE POTENTIALS 



Bp(p') 



B,{p) 



N(-p') 



N(-P) 



FIG. 1: The considered interaction in deriving the potentials. B\ (-B2) niay be Ac, Ec, or E*. The dashed line denotes the 
exchanged meson tt, a, p, or oi. 

The one-boson-exchange diagram in Fig. [T]is considered, where we use Bi and B2 to denote the charmed baryons. 
When deriving the potentials, we use the heavy quark limit for the charmed baryons, i.e. ignoring the l/A/gg 
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corrections but keeping up to corrections for the nucleon. At each interacting vertex, we introduce a cutoff A 

through the monopole type form factor 



A2 - g2 ' 



(12) 



where m is the mass of the exchanged meson and q is its 4-momentum. Because the meson exchange occurs between 
Ught quarks, we use the same cutoff for BiB2Tt and NNtt vertices. The cutoffs are taken to be around 1 GeV, while 
they may be different for the scalar and the vector meson exchanges. 

For convenience, we define some functions for the final potentials after the Fourier transformation. 
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(13) 



where k =^ p + p' with the notation in Fig. [l] Oi is the (transition) spin operator of the charmed baryon, O2 is the 

spin operator of the nucleon, and Oten is the tensor operator defined by Oten = '^'^'-'^''^^^'-^^'^^ — (^Q^ . 02)- The other 
definitions are 
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When getting these formulas, we have adopted the approximation q° = 0. The lowest excited state of AciV is the 
three-body KcNtt state. But in this study we do not consider the coupling of this channel because we need to excite 
the system to L = 1 in order to get the correct parity. Such excitation will hinder the channel coupling. 

In the multichannel case, we have transition potentials, e.g. V{KcN — > Yi*N), and we need to define the transition 
spin. Here we define u^^ = S'f <&, where 77'' is the Rarita-Schwinger field, is the transition spin matrix (2 x 4), and 
the spin wave functions of S* are 



$(3/2) = (1,0,0,0)^, $(1/2)^(0,1,0,0)^, $(-1/2)^(0,0,1,0)^, $(-3/2)^(0,0,0,1)^ 



(15) 



We also need the spin operator for S*, which is defined through (Jrs = ~Sl^aS^, Srs = f^rs- One may check 

Q2 — 3/3 I -1 \ 15 
^rs 2^2 ^ ' 4 ■ 

Now we may write down the general form of the potentials. 
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where i and j are labels of the channels, C{i,j) and C{i) are channel-dependent coefficients, the operator (T2 or O2 
is the Pauli matrix of the nucleon spin, and Oi is the corresponding spin matrix for the charmed baryon ai {(Jrs) or 
the transition spin St or S}. The uj exchange potential is similarly defined as the p potential. 

It should be noted that we omit the contact (^-functional) part of the potential, which appears in Y3 in our model. 
The main reason is that the interacting two baryons should get small contributions from such terms when they form 
a molecule-type bound state, in which two baryons are well separated. The difference between the defined functions 
Ys and Yi is the contact part. So the Ya does not appear in our model potentials. 

For the calculation with the potentials, one may use any consistent phase convention of the coefficients C{i,j). 
The correct result is ensured by Eq. ([2]). One should also note the possible convention problem when calculating 
the matrix elements {Oi ■ O2), (Oten), and {L ■ O) with the Wigner-Eckart theorem. The conventions must also be 
consistent. 



IV. THE COUPLING CONSTANTS 



The next task is to determine the coupling constants in the Lagrangians. We use several methods to constrain the 
values: (1) the strong decay of the baryons, (2) quark model estimation, (3) vector meson dominance assumption, 
(4) chiral multiplet assumption, and (5) QCD sum rule calculations. In Ref. [i^, the coefficients of BqBqtt, BqB^tt, 
B^Bq-k, BqB^tt, and BqBqTt coupling terms are called gi, 52, 93, 54, and 175, respectively. In the following parts, we 
also use these notations. 

The available strong decays are only Sc — > AcTt and S* — > AcTt [46| . The coupling constants 52 and 94 may be 
derived from the decay widths 

r(E,^A^) = ^^|p3|, r(S!^A^) = ^l^^|p3|, (17) 

where = 94 (180) MeV for Sc — > AcTt (E* — > AcTt) is the cm. momentum of the pion. In the numerical evaluation, 
we will use the averaged values derived from different decay modes: 1^2 1 = 0.598 and I54I = 0.999. 

All the coupling constants may be estimated with the quark model. The basic procedure is to calculate the average 
of the coupling vertex twice, at hadron level and at quark level, and to equate them. The effective Lagrangian for the 
light constitute u, d quark fields, -0 is 

^1 = -^^ri'd^i^'r'l^ - 9l^^^ - glh^'ip't.r' + ^t^)^ " Pp^a^^'' d^{f^y + c.,)^. (18) 

We assume that the heavy quark does not couple to the light mesons. For the pion coupling with heavy baryons, we 
just need to estimate two independent coupling constants, e.g. gi and 92, and relate others to them with the heavy 
quark symmetry relations 53 = -^51, <74 = — •\/3ff2, and 55 ~ —§51 [43| . The obtained values are gi = ^g\ and 
f/2 ~ ^\I\9'a- For the scalar meson coupling, we have £b = and £5 = 2g'^. For the vector meson coupling, one 

gets {Psgv) = -2g% [Psgv) = ^g% (Xsgv) = ^ + 4/^^ and (Xigv) = +4/p^), where m, is the quark 

mass. We have neglected the I/Mb^ ^ terms in deriving the above relations. For the nucleon-nucleon sector, we have 
9A = K = ig%, hy = g^, and hr - - |) + f 

With g\ = 0.75 in the chiral quark model [43], one obtains gi — 1.00, 52 = —0.61, (73 = 0.87, .94 = 1.06, g^ = —1.50, 
and gA = 1.25. For the coupling constant p*, the value in a ct model is 3.65 [i^. We have £3 = —3.65, £3 = 7.30, 
and ha = 10.95. If we use the value g^ = 2.621 in a chiral quark model [4^, we have £b = —2.621, is = 5.242, and 
h„ = 7.863. For the vector meson-baryon coupling constants, we use the quark mass niq = 313 MeV |4^ for the 
discussions. In the Nijmegen model D [H^l, the authors used hy = 2.1 and Ht = 9.1 GeV~^, which correspond to 
59 = 2.1 and = 2.8 GeV"! and gives {Pegv) = -4.2, (Psgv) = 8.4, {Xsgv) = 24.6 GeV"! and (Xigv) = -8.7 
GeV^^. In Ref. [5l[, the values = 3.0 and — 0.0 were used, which correspond to hy — 3.0 and hx ~ 6.4 GeV~^ 
and give (/Js^y) = -6.0, (Psgv) = 12.0, (Xsgy) = 19.2 GeV"!, (Xigy) = -6.8 GeV"!. 

In the heavy meson case, the relevant vector meson coupling constants have been estimated with VMD j52| . One 
can also use this approach to estimate the coupling constants /3b and f3s here. For the heavy quark hadrons, VMD 
assumes that the coupling of the electromagnetic current to the light quarks is dominated by the vector mesons. 
One may find the electromagnetic interactions of heavy hadrons in Ref. [53| . To avoid the confusion due to the 
existence of the heavy quark spectator, one uses only the isovector p meson dominance. The obtained values are 



6 



[Psgv] ~ — ^p'' = —5.04 and {(isgv) ~ — 10.08, where the decay constant fp = 216 MeV. Comparing these 

numbers with the quark model estimation, one finds that they are roughly consistent. 

In Ref. [U, the chiral partner of the ground charmed meson multiplet was studied to interpret Dsj(2317). In that 
method, one may obtain an estimation of the coupling constant between the scalar meson and the charmed mesons. 
The authors also explored the doubly heavy baryons. Recall that the parity partner of the nucleon has also been 
explored [s^. One may extend the previous studies to the charmed baryon case. Here, we derive the scalar meson 
coupling constants Ib and with the chiral multiplet assumption. 

For the antitriplet B^, its interpolating current is 

Bg ~ (<Z^C7^g)Q = {qlCqL)Q - {qlCqn)Q, (19) 

which transforms as (3,1) ® (1,3) under SU{i)L x SU{3)r. Here, we have used j^qL = qL and 7,5q_R = —qn- Its 
parity partner is 

^ {q^Cq)Q^{qlCqL)Q + {qlCqR)Q. (20) 

Therefore, one may define the left field iJg^ = -^[B^ + B^) and the right field = -^[B^ — Bg), which transform 
as _Bg^ — > LB^j^L^ and Bg^ — > RB^j^EF . Then the lowest order chirally invariant Lagrangian is: 

C?LH = \tr[B-^L{^l^~m-^)B-^:} + Ur[B-^^{zl^~m-^^^^^ (21) 

where S — ^cr^, E-^ is its transposition, and E transforms as E LUR^ . a can actually be the scalar nonet jH^j. In 
the chiral symmetric phase, the left and right fields have the degenerate mass mg. To get the linear representation, 
one redefines the fields Bg and Bg, 

Bg^ = -i=e(Bg + Bg)e^, B-^j,^ ^ {B-^ - B-^)C , (22) 

where the physical heavy field Bg or _Bg transforms as _B — > UBU^. With the fact that (a) « /^r in the chiral broken 
phase, we get 

^Ih = \tr[B-^iIpB-^ - rngSgSg] + Ur[B-^ipB-^ - rn-J-^B-^] + trlB^j^A^Bs + B-^^^Af^B-^] 

+9^tr{B-,B-, - BM + ^^tr{B-^aj,hyh - B^^phyBs} + ... (23) 
Now the mass difference AMg between i?g and Bg appears. Comparing with Eq. (jS]), one has 

AM, = G.fl i^^-^ = -^. (24) 

To estimate the numerical value of £b, we adopt the masses of the Ac states. There are three P-wave fields: Ad, Aco 
and Aci [EHl. The parity partner of the S-wave Ac should be Aco. Unfortunately, it has not been reported. We may 
use the A baryons to estimate the mass splitting AMg « A(1670) - A(1405) + Ac(2595)+ - A+ « 573 MeV, which 
gives £b = —3.1. This value and the sign are consistent with the quark model. 

The chiral multiplet (| , | ) of 5*^ (i''', |^) has not been identified. But we can still give a rough estimate 
of the coupling constant £s- Here, we discuss only the J = ^ parity partners Bq and Bq. If we use the interpolating 
current Bq ^ iq^Cjij,q)j^j^Q [13, [El], its partner's current is Bq ~ {q^Cjfj,j^q)^f^^Q, where = 7^ — w^^. We 
may define 

Bqrl - (qlCj^^qRHl'Q. Bqlr ~ {qlCj^^qLHl'Q. (25) 

In fact, {qRCjf.qL)'^ = {qlll^^^lR) = (ql'^^tJ^QR)^ so that Bqrl = Bqlr ^ (3,3), 

Bq ~ (q'^Cj^qhh'Q = {qRC^.qLHl'Q + (qlCl.qRhh'Q + ~ (3, 3), (26) 
Bq ^ (g^C7^7''Z)7tVQ = (9flC7^9L)7tV0 - (grC7^to)7r7'0 + ~ 0. (27) 
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This means the defined current Bq with the transformation (3,3) does not have a parity partner. However, one may 
use the following operators for the argument: 

Be ~ {q^Ca^,q)a^''Q ^ {qlCa^,qLh5CT''''Q - iqlCa^,qRh5'T^''Q, (28) 
Be - {q'^Ca^,q)^^a^''Q = {qlCa^^qLHa^^Q + {qHCa^,qRh5<y''''Q, (29) 

which transform as (6, 1) (1, 6). Then we may identify: 

They transform as BeL — > LBehL^ , BeR — > RBeR^F ■ With the similar procedure, we finally get is = ^J^- The 

parity partner of E is denoted as Eci [56| . It has not been measured, either. However, if we assume the mass 
difference comes mainly from the excitation of the light diquark, one may estimate AA/g w AA/g sa 573 MeV, and 
thus Is w 2£b « 6.2. 

Now we turn to some QCD sum rule calculations about the coupling constants. For the NNa coupling constant, 
there is a calculation ~ 14.6 [5^, which is a little larger than the estimation in the quark model. From the 
pion couplings S^S^tt and S*Ac7r in Ref. [e^, one gets 53 = 0.90 ± 0.17 ± 0.17, 54 = 0.94 ± 0.06 ± 0.20, and 
therefore gi « 1.04. From the values gsjs.Ti- = 4.2 ± 0.5 GeV"^ and ssja.tt = 7.8 ± 1.0 GeV~^ in Ref. [6l|, one has 
(73 = 0.78 ± 0.09 and 34 = 1.02 ± 0.13. These numbers are all consistent with the quark model results. In Ref. [6^ . 
however, the authors use a different definition of the coupling constant gs^Scir = —8.0 ± 1.7. With the help of the 
Goldberger-Treiman relation (GTR) [6^, one obtains a value gi ~ 0.3, which is much smaller than the quark model 
estimation. Ref. [5^ presents a calculation of the couplings of the light vector mesons and heavy baryons. To compare 
their results with our quark model estimation, we again use GTR and get (/3_b5\/) ~ —43.4 for gt^^t^^u = —1.85 ±0.06, 
{Psgv) ~ 123.6 for g^.^.^ = 2.65 ± 0.20, and {Psgv) ~ 78.0 for g^,^,^ = 1.51 ± 0.18. These results are much 
larger than the quark model estimation. But, for another coupling constant, we do not need GTR and we have 
(Xsgv) ~ 21.0 GeV-i for gl\^,p = 2.27 ± 0.20 and (Xsgv) « 13.5 GeV"! for g^\^,^ = 1.32 ± 0.18, which are not 
far from the quark model estimation. It is an interesting problem for the QCD sum rule calculation that we get 
inconsistent results when comparing different definitions of the coupling constants with the help of GTR in the heavy 
quark baryon case. Because of this reason, we do not use these values. By the way, one may conclude that the hidden 
local symmetry is supported by taking a glimpse at the small coupling constants g^,^,p and g^,^*i^ in Ref. [58| . 
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TABLE II: The coupling constants in different methods. For the quark model estimation, we present the results with g\ = 0.75, 
Qa ~ 3.65 [4^ . Qp — 3.0, and = 0.0 [5l[. We have used the phase convention of g2,4 from the decay being consistent with 
the quark model. 



For comparison, we summarize in Table HIl the values of the coupling constants given by different methods. There 
we have taken our convention for the phases. In our numerical evaluation, we will choose the following values and 
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relations, 

0.598, 54 - 0.999, 91 = ^gi, .93 = ^154, 35 = -V2r74, 5a = 1.25, 

= -3.1, is = -2^s, K = 10.95, 
-6.0, (fisgv) ^ -'2{I3b9v), (Asgy) = 19.2 GeV-\ [Xigv] = -{^s9v)/ ^/^, 

/iy = 3.0, /iT = 6.4GcV"^ (31) 

V. = 0+ CASE 

We adopt the variational method [i^l to solve the bound state problem. For the hadron masses, we use = 137.27 
MeV, mp = 775.49 MeV, = 782.65 McV, tun = 938.92 McV, mA, = 2286.46 MeV, ms, = 2453.56 MeV, and 
ms* = 2517.97 MeV [i^. The scalar a meson is a broad state with strong coupling to tttt S-wave scattering states 
and its mass has not been accurately given. The recent investigations indicate that the pole mass is around 400 ~ 600 
MeV [65l-[67j. In the following calculation, we just adopt a large value m„ = 600 MeV. If we take a smaller value, 
then it will enhance attraction and will result in deeper bound states. Another parameter we cannot determine is the 
cutoff. In fact, the results are sensitive to it. We treat it as a free parameter and discuss the effects on the binding 
energy. In this section, we investigate the spin-singlet case. 

It is useful to take a look at the potentials. We plot the potentials of different channels in Fig. [5] with the cutoff 
= Ao- = Avec = 1 GeV. The pion potentials in the channels 2 and 3 are both repulsive while the transition potentials 
are strong. The central force of the p potential dominates the transition process KcN{^So) 'T.cNQ-Sq) while the 
tensor force of the tt potential dominates the processes KcNQ-So) o Yi*N{^Dq) and TicN{^Sa) <-> S*iV(^£)o)- We 
study first whether only one-pion-exchange interaction would lead to a bound state and then investigate the scalar and 
vector meson contributions. We call the corresponding potentials OPEP (one-pion-exchange potential) and OBEP 
(one-boson-exchange potential), respectively. Certainly the portions of the contributions from different mesons change 
if one uses other cutoffs A^ 7^ A^ 7^ Avoc- We will discuss the effects in the following calculation. 

A. OPEP model 

Let us first consider the case without channel coupling. For the channel AciV, the direct one-pion-exchange is 
forbidden. For the channels ScTV and S*7V, the potentials are both repulsive, therefore one cannot get a bound state. 



A. (GeV) 


1.2 1.3 


1.4 


1.5 1.6 1.7 


B.E.{J = 0) (MeV) 


0.64 6.16 


18.51 


38.88 68.29 107.64 


VW) (fm) 


5.2 1.9 


1.2 


0.9 0.8 0.7 


Prob. (%) 


(98.2/0.6/1.2) (94.0/2.3/3.7) ( 


89.3/4.6/6.1) ( 


84.5/7.2/8.3) (80.1/9.8/10.1) (76.1/12.2/11.7) 



TABLE III: Binding solutions for the = 0^ case with channel couphng in the OPEP model. The binding energies (B.E.) 
are given relative to the Ac A'" threshold. The probabilities correspond to AcA*'(^5'o), T,cN(^ So), and E*A''(''_Do), respectively. 

For the case with channel coupling, we have one free parameter A^ . It is interesting that a bound state is obtained 
within the reasonable range of the cutoff. We present the binding energies, the root-mean-square (RMS) radius, and 
the probabilities of each channel in Table IIIII The binding energy is given with relative to the Ac V threshold. As a 
deuteronlike molecule, the two baryons should not be very close. Thus wc only list the results with the RMS radius 
larger than 0.7 fm. The values indicate that the reasonable binding energy should be no more than tens of MeV 
although it depends on the poorly known cutoff. As an example, we show the wave functions of different channels 
for A^ = 1.3 GeV in Fig. [3] (a). That the probability of the channel E*iV(^Do) is larger than that of the channel 
TicN{}Sq) indicates the importance of the tensor force in the model. As a check, we have calculated the two channel 
case: AcA^(^S'o) and ScA^(^S'o) and we do not find any binding solutions. 

We have omitted the (5-functional part in our potentials. Once that part is included, deeper molecular bound states 
arc obtained. In that case, the S-wave TjcN channel is dominant over to D-wave Tf^N channel. For example, if we 
use Att = 0.8 GeV, a 3-channel calculation gives the binding energy B.E. = 19.18 MeV and a 2-channel calculation 
(without the _D-wave channel) gives B.E. = 17.85 McV. However, the binding energy is more sensitive to the cutoff 
parameter. We get B.E. = 111.53 McV in the 3-channel calculation with a little larger cutoff A^r = 0.9 GeV. 



52 = - 
{Psgv) = - 
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FIG. 2: The potentials of difTerent channels for the = 0"*" case with = An 



Av 



1 GeV. 



B. OBEP model 



After we consider the contributions from the scalar meson a and the vector mesons p and oj, the binding energies 
will change accordingly. However, wc have some freedoms to choose: the three cutoff parameters A^, Aq-, and Avoc- 
In the nucleon-nucleon case, one may determine the unknown model parameters by fitting the abundant experimental 
data. Since there is no available data, here we discuss the results with two approaches for the parametrization of the 
cutoffs. 
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1. Common cutoff 



Acom (GcV) 


0.9 1.0 1.1 1.2 1.3 1.4 


A,iV(^5'o) [B.E.., 
T,,NCSo) [B.E., ypy] 


[1.24, 3.8] [11.09, 1.5] [27.07, 1.1] [46.66, 0.9] [68.45, 0.8] [91.58, 0.7] 
X [2.22, 2.8] [14.22, 1.3] [31.56, 1.0] [52.08, 0.8] [74.62, 0.7] 



TABLE IV: Binding solutions for the individual channels in the = 0^ and common-cutoff case in the OBEP model. The 
binding energies (B.E.) are given relative to their own thresholds, "x" indicates that there is no binding solution. The units 
for the binding energy and RMS radius are MeV and fm, respectively. 



Ae 



(GeV) 



B.E.{J = 0) (MeV) 

ypy (fm) 

Prob. (%) 



0.8 



0.9 



1.0 



1.1 



0.12 13.60 52.50 123.14 

11.2 1.5 0.9 0.7 

(99.7/0.1/0.2) (96.0/2.0/2.0) (87.3/9.2/3.5) (75.8/19.7/4.5) 



TABLE V: Binding solutions for the = 0+ and common-cutoff case with channel coupling in the OBEP model. The binding 
energies (B.E.) are given relative to the AcA^ threshold. The probabilities correspond to AcA'^(^S'o), TicNC^ So), and T,'^N{^ Do), 
respectively. 



The simplest assumption is that we use the same cutoff for different mesons A^r = = Aypc = ^com- This cutoff 
should be at least larger than the exchanged meson masses. As in the former subsection, we first consider the case 
without coupled channel effects. In Table IIVI we show the binding energy and the corresponding RMS radius for 
individual channels relative to each threshold, where solutions are found only in the S-wavc channels. Numerically, 
one may obtain binding solutions for the D-wave channel with a larger cutoff, but the solutions arc not reasonable. 
For example, wc get B.E. = 1.99 McV and trms = 0.6 fm with Acom = 1-78 GeV. The radius is so small for a 
shallow D- wave bound state. Because of the attractive a meson, AcN can be bound now. The contribution from 
the repulsive lj is not large in this parametrization of cutoffs. For the 'ScN , the cancellation between the p and uj 
contributions is large and the attraction comes mainly from the scalar potential. 

Now we consider the coupled channel effects. Table |V] lists the results. The binding energy and the radius are 
more sensitive to the change of the cutoff than the OPEP case. For the same binding energy as the OPEP case, the 
necessary cutoff in OBEP is smaller. These features are attributed to the more attractive potentials. In order to have 
the molecular condition, trms > 0.7 fm, satisfied, the binding energy should be no more than tens of MeV, as in the 
OPEP model calculation. As an example, we plot the wave functions of different channels with Acom = 0.9 GeV in 
Fig. Kb). 

Let us go back to the single channel AcA'^ and the potentials. Because of the isospin conservation, direct tt and p 
exchanges are forbidden. The spin-dependent a; exchange interaction also vanishes. The resulting spin- independent 
repulsive potential has been shown in Fig. O If one includes the (5-functional terms in the model potentials, the u! 
exchange interaction is strongly attractive at short distance, which results in much deeper bound states and enhances 
the sensitivity of the results to the cutoff parameter. In addition, the coupled channel calculation with the (5-fmictional 
terms results in unreasonable molecular bound state solutions. Thus we omit the J-functional terms in the present 
model construction. 



2. Scaled cutoffs 

Another possible choice of the cutoff parameters is to choose different values for the pseudoscalar, scalar, and vector 
mesons. To reduce the number of parameters, we adopt a parametrization used in Ref. (68j . Aox = JtIcx + ctAgcD- 
Here TOox is the mass of the exchanged meson, Aqcd = 220 McV is the scale of QCD, and a is a dimcnsionless 
parameter whose value is not very far from 1. We choose various values of a in the present investigation. 

Table IVTl shows our results for the case without channel coupling. For the D- wave channel E* A^, one gets a solution 
with the binding energy 0.92 MeV and the RMS radius 0.7 fm by using a = 3.54. A larger value a = 3.6 gives a 
solution with the binding energy 21.28 MeV and the radius 0.6 fm. Therefore we concern only S-wave states in that 
table, as in the case of common cutoff. By comparing these values with those in Table llVi one finds that these two 
parameterizations give similar results, especially for the l^cN channel. 

For the case with channel coupling, we give the numerical results in Table IVTIl From this table and Tablc|Vl although 
the binding energies and the radii in these two parametrization methods are consistent and somehow equivalent, while 
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A^NCSo) [B.E.,y/{r 



1.5 



2.0 



2.5 



3.0 



3.5 



4.0 



[0.12, 11.6] [6.54, 1.9] [20.30, 1.2] [38.86, 0.9] [60.56, 0.8] [84.29, 0.7] 
X [2.33, 2.7] [14.07, 1.3] [31.64, 1.0] [52.96, 0.8] [74.75, 0.7] 



TABLE VI: Binding solutions for the individual channels in the = 0"'" and scaled-cutoff case in the OBEP model. The 
binding energies (B.E.) are given relative to their own thresholds, "x" indicates that there is no binding solution. The units 
for the binding energy and RMS radius are MeV and fm, respectively. 



B.E.{J = 0) (MeV) 

\fW) (fm) 
Prob. (%) 



1.2 



1.4 



1.6 



1.8 



2.0 



0.11 5.26 
11.7 2.0 
(99.6/0.4/0.0) (95.8/4.1/0.1) 



19.37 43.05 75.65 

1.2 0.9 0.7 

'.7/10.1/0.2) (76.9/22.5/0.6) (77.0/22.5/0.5) 



TABLE VIL Binding solutions for the — O""" and scaled-cutoff case with channel coupling in the OBEP model. The binding 
energies (B.E.) are given relative to the Ac A'' threshold. The probabilities correspond to AcA(^S'o), ScA(^5'o), and E*A''('''_Do), 
respectively. 



the probabilities of individual channels are not. The probability depends on the contributions from different exchanged 
mesons. It is helpful to understand this feature from the wave functions and the potentials. We show the wave functions 
of different channels with a = 1.4 in Fig. [3] (c). There is a node in the wave function of the third channel 'E*N{^Do). 
The node does not disappear for the other a. Actually, such a behavior results from the transition potential between 
the first channel and the third channel as well as that between the second channel and the third channel. These 
potentials have different signs in the long-range part and the short-range part. As an example, we plot the potentials 
of this parametrization in Fig. |3]with a — 1.4. 



VI. = 1+ CASE 



The spin-triplet case is a little complicated because there are four more channels than the spin-singlet case. We 
carry out similar analysis in this section. Before the numerical evaluation, we plot the potentials with the cutoffs 
Att = Act = Avoc = 1 GeV, which are shown in Fig. [S] The potentials for the channel AcN{^Si) O AcN{^Si) and 
the channel AcNC-Sq) O AcN(^So) are identical. There is no transition potential AcA^(^S'i) O AcN{^Di), which is 
omitted in Fig. [51 



A. OPEP model 



Now we discuss possibilities of bound states without channel coupling among AcN, T,cN, and T^^N in the OPEP 
model. For the AcN channel, there is no binding solution as before. For the 'ScN and T,*N channels, S- and D- 
wave mixing is allowed. The former (latter) is a two-channel (three-channel) problem: EcA^('^S'i) and T^cN{'^Di) 
(S*Ar(35i), E*N{^Di), and E*N{^Di)). Numerical results are given in TableEnH 



A, (GeV) 




1.0 1.1 1.2 1.3 


1.4 


1.5 


1.6 1.7 1.8 


B.E.{J = 1) (MeV) 




0.08 1.22 4.09 9.27 


17.34 


28.93 


44.69 65.32 91.52 


ypy (fm) 




13.8 4.0 2.4 1.7 


1.4 


1.1 


1.0 0.8 0.7 


Prob. (%) 


(9J 


S.7/1.3) (95.8/4.2) (93.5/6.5) (91.7/8.3) (90.2/9. 


S) (89.0/11.0) 


(88.0/12.0) (87.2/12.8) (86.4/13.6) 


A, (GeV) 


1.6 1.7 1.8 




1.9 


2.0 2.1 


B.E.(J = 1) (MeV) 


0.47 3.83 10.27 




20.08 


33.66 51.50 


S*A (fm) 




5.6 2.2 1.5 




1.1 


0.9 0.8 


Prob. (%) 




(96.3/0.1/3.6) (92.9/0.1/7.0) (91.0/0.2/8. 


B) (89 


7/0.2/10.1) (8 


5.8/0.1/11.1) (87.9/0.1/12.0) 


D-wave prob. 




3.7% 7.1% 9.0% 




10.3% 


11.2% 12.1% 



TABLE VIII: Binding solutions for the individual channels in the = l""" case in the OPEP model. The binding energies 
(B.E.) are given relative to their own thresholds. The probabilities correspond to "^Si and ^Di for the EcA system, and "^Si, 
^D\, and ^Di for the E*A'' system, respectively. We also present the total D-wave probability for E*A. 



After including the coupled channel effects, we obtain the results presented in Table HXl Our calculation indicates 
again that the coupled channel effects are important in the OPEP model. One finds that the binding energies are 
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(GeV) 


1.2 1.3 1.4 1.5 1.6 


B.E.{J = 1) (MeV) 
(fm) 
Prob. (%) 

D-waye prob. 


1.17 7.52 20.64 41.86 72.39 
3.9 1.8 1.2 0.9 0.8 
(97.6/0.2/0.4/ (93.7/0.6/1.3/ (89.3/1.4/2.5/ (84.6/2.7/3.7/ (79.6/4.6/4.9/ 
0.0/1.0/0.1/0.7) 0.0/2.5/0.2/1.7) 0.0/3.9/0.3/2.6) 0.1/5.3/0.4/3.2) 0.1/6.8/0.4/3.6) 
1.8% 4.4% 6.8% 9.0% 10.1% 



TABLE IX: Binding solutions for the = 1^ case with channel coupling in the OPEP model. The binding energies (B.E.) are 
given relative to the A^N threshold. The probabilities correspond to AcAf(^Si), EciV(^5'i), 'S*N{^Si), A^Ni^Di), EcAf(^Di), 
E*A''(^_Di), and T,lN{^Di), respectively. We also present the total D-wave probability. 



slightly different between the singlet and triplet cases. The 4-th channel, AcN{'^ Di), has small contribution. Actually, 
if one ignores this channel, the resultant binding energy changes little. We show the wave functions of different 
channels with the cutoff A^ — 1.3 GeV in Fig. [3] (d). 



B. OBEP model 



In the case of the OBEP model, we again explore two cases for the parametrization of the cutoff parameters: (1) 
common cutoffs A^r = = Avoc = Acom, and (2) scaled cutoffs Acx — ruex + cxAqcd- 



1. Common cutoff 



We first discuss the case without the channel coupling. For the AcA^ channel, the scalar meson and the vector 
meson exchanges do not lead to S- and D-wavc mixing and thus one again gets the binding solutions in Table IIVI 
i.e. the binding solution for the single channel AcA^ is spin independent. For the TicN and Yi*^N channels, we present 
their results in Table [Xj 



Acom (GeV) 


0.8 0.9 1.0 






B.E.{J = 1) (MeV) 


2.84 26.20 78.51 




ScA (fm) 




2.8 1.2 0.8 






Prob. (%) 




(96.4/3.6) (95.0/5.0) (95.0/5.0) 




Acom (GeV) 


0.9 


1.1 1.3 1.5 


1.7 


B.E.{J = 1) (MeV) 


0.73 


13.10 27.04 42.89 


77.36 


S^A y(r2) (fm) 


4.7 


1.4 1.1 0.9 


0.8 


Prob. (%) 


(97.7/0.1/2.2) (94.4/0.1/5.5) (91.3/0.1/8.6) (85.1/0.0/14.9) (63.5/0.5/36.0) 


D-wave prob. 


2.3% 


5.6% 8.7% 14.9% 


36.5% 



TABLE X: Binding solutions for the individual channels in the — 1"*" and common-cutoff case in the OBEP model. The 
binding energies (B.E.) are given relative to their own thresholds. The probabilities correspond to ^Si and ^Di for the EcA 
system, and "^Si, ^Di, and ''Di for the E*A system, respectively. We also present the total D-wave probability for E*A. 



After considering the coupled channel effects, we get the results in Table IXll which are consistent with the OPEP 
model calculation. An example of the resulting wave functions with Acom = 0.9 GeV is given in Fig. [3] (e). 



Ac. 



(GeV) 



B.E.{J = 1) (MeV) 
(fm) 
Prob. (%) 

D-wave prob. 



0.8 



0.9 



1.0 



1.1 



0.22 
8.7 

(99.6/0.0/0.1/ 
0.0/0.2/0.0/0.1) 
0.3% 



13.49 
1.5 

(96.6/0.3/0.9/ 
0.0/1.2/0.1/0.9) 
2.2% 



47.50 
0.9 

(90.6/2.2/3.7/ 
0.0/2.1/0.1/1.3) 
3.5% 



106.16 
0.7 

(79.1/9.9/6.4/ 
0.1/3.2/0.1/1.2) 
4.6% 



TABLE XL Binding solutions for the = 1+ and common-cutoff case with channel coupling in the OBEP model. The binding 
energies (B.E.) are given relative to the AcA threshold. The probabilities correspond to AcA(^S'i), EcA*'(^5'i), E*iV(^5'i), 
AcA(^Di), EcA(^Di), E*A(^Di), and E*A('^Di), respectively. We also present the total D-wave probability. 
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2. Scaled cutoffs 



a 


0.9 1.1 1.3 1.5 1.7 1.9 


B.E.{J = 1) (MeV) 

E,iv ypy (fm) 

Prob. (%) 


0.39 6.37 19.52 39.26 64.90 95.77 
6.6 2.0 1.3 1.0 0.8 0.7 
(99.8/0.2) (99.7/0.3) (99.8/0.2) (99.8/0.2) (99.9/0.1) (99.9/0.1) 


a 


3.2 3.3 3.4 3.5 3.6 3.7 


B.E.{J = 1) (MeV) 
(fm) 
Prob. (%) 
D-wave prob. 


1.39 3.13 6.39 12.41 22.74 38.38 
3.4 2.4 1.7 1.3 0.9 0.7 
(98.0/0.1/1.9) (96.4/0.1/3.5) (92.8/0.2/7.0) (85.3/0.2/14.5) (73.7/0.2/26.1) (61.1/0.3/38.7) 
2.0% 3.6% 7.2% 14.7% 26.3% 39.0% 



TABLE XII: Binding solutions for the individual channels in the — 1^ and scaled-cutoff case in the OBEP model. The 
binding energies (B.E.) are given relative to their own thresholds. The probabilities correspond to '^Si and ^Di for the J^cN 
system, and "^Si, ^-Di, and ^Di for the E*A'^ system, respectively. We also present the total D-wave probability for T,*N. 



We present the results without channel coupling in Table IXIII and those with channel coupling in Tabic IXIIII If one 
compares Table IXlIIl with Table IxH it is easy to see that the two parameterizations of the cutoffs may give consistent 
binding energies and radii. The wave functions for a = 1.5 in the coupled channel case are given in Fig. [3](f). As in 
the corresponding spin-singlet case, the wave functions also have nodes, which are due to the transition potentials. 
In fact, the nodes have appeared in the D wave of the S*Af channel (without the channel coupling to the AcN and 
S,iV). ' 



a 


1.3 1.5 1.7 1.9 2.1 


B.E.{J = 1) (MeV) 

VW) (fm) 
Prob. (%) 

-D-wave prob. 


0.66 6.63 19.82 40.98 70.27 
5.1 1.9 1.2 0.9 0.7 
(99.3/0.2/0.5/ (96.7/1.0/2.2/ (92.5/2.8/4.5/ (86.9/5.8/6.8/ (80.5/10.2/8.6/ 
0.0/0.0/0.0/0.0) 0.0/0.0/0.0/0.1) 0.0/0.1/0.0/0.1) 0.0/0.2/0.0/0.3) 0.0/0.3/0.0/0.4) 
0.0 0.1% 0.2% 0.5% 0.7% 



TABLE XIII: Binding solutions for the — 1^ and scaled-cutoff case with channel coupling in the OBEP model. The binding 
energies (B.E.) are given relative to the AcA'^ threshold. The probabilities correspond to AcA'^(^S'i), EcA''(''Si), E*A''('^5'i), 
K^Nl^Di), T.cN{^Di), E*A'(^Di), and E*A'(^Di), respectively. We also present the total D-wave probability. 



VII. SUMMARY AND CONCLUSIONS 

In studying the possible molecular bound states containing AcN, we have constructed a potential model based on 
the effective Lagrangian reflecting the heavy quark symmetry, chiral symmetry, and hidden local symmetry. Solving 
the Schrodinger equation for the coupled AcA^ — Y,cN — Y,*N systems for J = and J = 1 states, we obtain molecular 
bound states in both channels with appropriate cutoffs. By analyzing all the binding energies and the corresponding 
RMS radii, one observes several features: (1). the spin-singlet state and the spin-triplet state have slightly different 
binding energies for a given radius and the reasonable binding energy is at most tens of MeV; (2). the OPEP model 
and the OBEP model are somehow equivalent in getting the consistent binding energy and the corresponding RMS 
radius; (3). the coupled channel effects are important both in the OPEP model and in the OBEP model. To see these 
features, we make a comparison, which is presented in Table IXIVl 







AcN (S-wave) A^A^ - E^A - E*A 


0+ 


OPEP (A) 
OBEP (A) 
OBEP (a) 


X [1.367 
[0.900: 1.24, 3.86] [0.900 
[1.533: 0.25, 8.13] [1.533 


13.60, 1.38] 
13.60, 1.46] 
13.57, 1.37] 


1+ 


OPEP (A) 
OBEP (A) 
OBEP (a) 


X [1.353 
[0.900: 1.24, 3.86] [0.900 
[1.618: 0.80, 4.72] [1.618 


13.54, 1.40] 
13.49, 1.47] 
13.47, 1.39] 



TABLE XIV: Comparison between different cases. The meaning of the numbers are [cutoff A in GeV or dimensionless a: 
binding energy in MeV, RMS radius in fm]. 
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32 B.E. 
-0.50 13.17, 11.36 
-0.55 13.38, 12.38 
-0.60 13.61, 13.54 
-0.65 13.83, 14.84 
-0.70 14.07, 16.29 


Qi B.E. 
0.7 6.93, 9.54 
0.8 8.77, 10.67 
0.9 10.99, 11.99 

1.0 13.62, 13.51 

1.1 16.68, 15.26 


£b B.E. 
-2.0 2.09, 4.15 
-2.5 6.24, 9.15 
-3.0 12.21, 15.75 
-3.5 19.75, 23.76 
-4.0 28.68, 33.03 


iPBgv) B.E. 
-3.0 14.75, 14.67 
-4.0 14.36, 14.27 
-5.0 13.98, 13.88 
-6.0 13.60, 13.49 
-7.0 13.22, 13.11 


(Xsgv) B.E. 
13.0 13.03, 13.21 
15.0 13.20, 13.29 
17.0 13.39, 13.39 
19.0 13.58, 13.48 
21.0 13.78, 13.59 


h„ B.E. 
8.0 4.04, 4.10 
10.0 10.00, 9.96 
12.0 18.08, 17.90 
14.0 27.98, 27.63 
16.0 39.46, 38.93 


hv B.E. 
2.0 13.92, 13.89 
2.5 13.76, 13.69 
3.0 13.60, 13.49 
3.5 13.44, 13.30 
4.0 13.28, 13.10 


hr B.E. 
5.0 13.48, 13.53 
6.0 13.56, 13.50 
7.0 13.66, 13.48 
8.0 13.77, 13.47 
9.0 13.89, 13.47 



TABLE XV: J — and J = 1 binding energies (MeV) corresponding to the variation of only one coupling constant. The 
unchanged coupling constants are the same as in Eq. (31). The units of (Xsgv) and hr are GeV^^. The common cutoff 
Acorn ~ 0.9 GeV is used. 



The deuteron, the well-known pn bound state with the binding energy of 2.2 MeV, is also a spin-triplet state. The 
D-WAYC contribution is very important and the corresponding probability of the state at D-wave is around 4^6%. If 
the D-wavc probability of AcN were in this range, the binding energy would be around 10 MeV in the OPEP model 
from Table El 

For the A^A^ single channel, only cr and w exchanges are allowed. There are no tensor force nor spin-dependent 
parts in the potentials. As a result, the coupling between the ^Si and ^Di states vanishes and one has ViSo = V^Si 
for the potentials. That is, the possible = 0+ and = 1+ molecular states are degenerate in the heavy quark 
limit. After considering the contributions from the channels YlcN and 'E*N, the tensor force and the spin-dependent 
interactions enter. However, full coupled channel calculation does not give significantly different results for the two 
spin states. That is, the two spin states are still qualitatively degenerate. Maybe this feature indicates that the 
coupled channel effects do not change the degeneracy in the dominant channel. If this is the case, one expects that 
the states with different angular momenta of a two-body system like AcSc, AcA, AcS, or AcA would have similar 
binding solutions. They may be tested in the future investigations of bound state problems and scattering problems. 

The results are sensitive to the phenomenological cutoff parameter, which is a general feature for the molecule 
study in the meson exchange models at the hadron level. In effect, this parameter encodes the size effects of the 
hadrons and the information of the short-range interaction. This feature indicates that the binding energy is sensitive 
to the short-range interaction. At this moment, we cannot determine the values of the cutoff. So we treat it as a 
free parameter and discuss the results for a reasonable range of the cutoff. If the cutoff parameter around 1.2 1.4 
GeV in the OPEP model, 0.8 1.0 GeV in the OBEP model (common cutoff case), or 1.2 1.7 in the OBEP model 
(scaled cutoff case) is reasonable, then one obtains moleculelike bound states. Larger parameters may be beyond 
present models while smaller ones would result in unbound states. 

To see the sensitivity of the binding energy to the cutoff parameter A or a, we show the results in a diagrammatic 
form. In Fig. [HI we present the dependence of the binding energy on A^r for the = 0+ and = 1+ cases in the 
OPEP model. The dependence of the binding energy on the common cutoff Acom is illustrated in Fig. [T] That on the 
parameter a is shown in Fig. |S1 The binding energies in these figures arc not all given in the former tables. Here, we 
have included some additional binding energies whose corresponding RMS radii are smaller than 0.7 fm. From these 
figures, it is obvious that the coupling among the channels AcN, T,cN, and T,*N may result in bound states even if 
there are no binding solutions within individual channels. 

In the OPEP model, the attraction for the formation of the AcN bound states comes only through the T,cN and 
E*iV channels indirectly. So the coupled channel effects are essential. The same observation is made in the AcAc case 
In the OBEP model, the single channel AcN itself may have binding solutions. Then the effects of the channel 
coupling make the binding energy larger with a significant value. 

In the above discussions, we do not consider the uncertainty of the coupling constants in Eq. (31). To investigate 
the sensitivity of the results to them, we vary separately 52, 94, ^b, [Psgv), i/3s9v), {^s9v), ^o-, hy, and Ht and 
compare the binding energies. As a well-known value, gA is not varied. We still adopt the relations in Eqs. (31) in 
determining other coupling constants. Table IXVl illustrates the sensitivity when we use the common cutoff parameter 
Acom = 0.9 GeV. Larger values of the coupling constants (except {(3b9v) and hy) may result in deeper bound states. 
From the numerical results, one finds that the results are sensitive to (74, ^b, and h^- If one uses the set of coupling 
constants 32 — —0.5, 54 = 0.7, £b = —2.0, [Psgv) = —7.0, {\s9v) = 13.0, = 8.0, hy = 4.0, and hr = 5.0, one may 
also get binding solutions with a slightly larger cutoff Acom = 1.1 GeV. Since this set is not helpful to the formation 
of bound states, the main conclusion that the molecular bound states are possible does not change. 
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Theoretically, the formation of these = O"*", l"*" KcN bound states occurs after the production of Ac and through 
the coalescence mechanism, similar to that of X(3872) in the molecular picture [T^I- Once the molecules are formed, 
they would be rather stable because they do not decay through strong interactions. Possible places to search for them 
are GSI-FAIR, J-PARC, and RHIC. It is also possible to look for them at BELLE. 

In summary, we have obtained effective Lagrangians describing the charmed baryon interactions with light mesons 
according to the heavy quark symmetry, chiral symmetry, and hidden local symmetry. We estimate the coupling con- 
stants with various methods and get consistent results. Based on these interaction Lagrangians, one-boson exchange 
potentials for the A^A^ — Sc-^ — S*7V systems arc derived. With the cutoffs around 1 GeV, it is possible to get binding 
solutions by solving the bound state problem with the variational method. One finds that the channel coupling has 
very important effects for the possible molecular bound states. These effects do not significantly change the feature 
that the binding solutions are spin-independent. If these states exist, the reasonable binding energies should be at 
most tens of MeV. 
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Appendix A: Phase convention in the superfleld 



If one reduces Eq. (3.12) in Ref. to the relevant terms in Eq. ([3]) in this paper with the superfleld defined 
in Eq. ([T|), one gets 



33 = --^^,91, .95 = -^.91, 94 = V3Sg2, .96 = 0. 



(Al) 



To determine the phase 5, we resort to the quark model and the heavy quark symmetry. The spin wave functions of 
and E* (5*^ = 5) are 



s*(t) 

E(t) 



(11); 



(11); 



(10) t, 



(10) 



(A2) 



where e is an arbitrary relative phase, ; or f is the heavy quark spin, and (11) or (10) represents the diquark spin 
and its z component. Let af^ denotes the spin operator acting only on the heavy quark, and 

<E*(t) = aE*(t) + feE(t), 
aim) = cE*(t)+dI](t). (A3) 



One gets a= ^, b 



heavy quark transformation. In the static limit, the superfield S*^; {fi = 3) is 



4=e, and d — — i. Now we demand the superfield 5*^ does not change under the 



'71 



(7. 



(A4) 



where and St are transition spin operators of Rarita-Schwinger field and Dirac field, respectively, $ = (0, 1, 0, 0)^ 
is the wave function of E*, and x = (0, 1, 0, 0)"^ is the wave function of Ec. With (jf^S^ — S3, one may get {6e) = —1. 
A natural choice is e = 1 and then 6 = —1. Therefore, we have 



.93 = ^91, 95 = -fffi, 54 = -^/3g2, 9e = 0. 
The relations among the coupling constants have been presented in Ref. [i^ . 



(A5) 
(A6) 
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Channel 





0,0 1.0 2.0 3.0 4.0 0.0 1.0 2.0 3.0 4.0 

r [fm] r [fm] 

(e) (f) 

FIG. 3: The wave functions Ui (i=l,2, 7) of different channels. The first three diagrams correspond to the results for the 
spin-singlet state: (a) OPEP case with the cutoff = 1.3 GeV; (b) OBEP case with the cutoff = = Avoc = 0.9 GeV; 
(c) OBEP case with the parameter a = 1.4. The last three diagrams correspond to the results for the spin-triplet state: (d) 
OPEP case with the cutoff A^ = 1.3 GeV; (e) OBEP case with the cutoff A^ = A<^ = Avcc = 0.9 GeV; (f) OBEP case with the 
parameter a — 1.5. 



19 





120.0- 




60.0 


, — , 


0.0 


>- 






-60.0 








-120.0- 








-180.0 




-240.0- 




-300.0- 



0.0 



1.0 

r [fm] 
(11): AciV(i5o) o AciV(iSo) 



2.0 



> 
CD 



1170.0- 
585.0- 
0.0- 
-585.0- 
-1170.0- 
-1755.0- 
-2340.0- 
-2925.0- 



9 

u 

Tot 



1.0 

r [fm] 



(33): E*iV('^Do) o YT.Ni^D^) 



> 



25.0- 
0.0- 
-25.0 -, 
-50.0- 
-75.0- 
-100.0- 
-125.0- 



0.0 1.0 

r [fm] 
(13): AcAf('So) o E*Arf A)) 




1.0 

r [fm] 
(22): E,iV(^5o) O E,A/-(iSo) 



2.0 



0.0 

-45.0 

-90.0 

-135.0 

-180.0 

^ -225.0 

-270.0 

-315.0 

-360.0 

0.0 



Tot 



1.0 

r [fm] 
(12): A,iV(iS'o) « E,Ar(i5'o) 



40.0 
20.0 
0.0 
-20.0 
^ -40.0 
-60.0 
-80.0 



0.0 1.0 

r [fm] 
(23): EciV('So) « E:iV(5Do) 



2.0 



FIG. 4: The potentials of different channels for the = 0^ case with a = 1.4. 



95.0 
0.0 
-95.0 
-190.0 
-285.0 
-380.0 
-475.0 
-570.0 











Meson 

7T 


cr 

P 

CJ 

■ ■ ■ ■ Tot 



0.0 1.0 

r [fm] 

(22): EcAf(^S'i) ^ EciV(^S'i) 



2,0 




1.0 

r [fm] 

(44): AciV(^I>i) ^ AcAf(^Di) 



2.0 



100.0 
0.0 

-100.0 
-200.0 
-300.0 
-400.0 
-500.0 



/ ^ - _ 
K ^ ^ - ~ ~ ~ ^ 








Meson 

7T 


a 

P 

CJ 

■ ■ ■ ■ Tot 



0.0 1.0 

r [fm] 

(66): T,*N(^Di) « T.*N{^Di] 



2.0 




0.0 1.0 

r [fm] 

(33): E*iV(^S'i) o E*iV(^S'i) 



2.0 



> 



370.0 
185.0 
0.0 
-185.0 
-370.0 
-555.0 
-740.0 



-925.0- 























■ ' / 






/ 






.' / 






■ / 






- ■ / 


Meson 




/ 


7T 




/ 


a 




/ 


P 






u 






Tot 





0.0 1.0 

r [fm] 

(55): ^cN^Di) ^ E,Ar(^Di) 



2.0 



> 
CD 




-1325.0 



1.0 

r [fm] 

(77): T.*NCDi) « ElN^Di) 



2.0 



FIG. 5: The potentials of different channels for the — 1^ case with A.^ = A^ = Avcc ~ 1 GeV. 



21 




1.0 

r [fm] 

(12): A,iV(=*5'i) ^ EeiV(=*5'i) 



2,0 



160,0 
^ 120,0 
2, 80,0 
40.0 
0,0 





Meson 

TT 


a 

P 

— - - CJ 
■ ■ ■ ■ Tot 











0,0 1,0 

r [fm] 

(15): AcNi^Si) ^ EeAf(3£'i) 



2,0 



40,0 
0,0 

r- -40,0 
> 

-80,0 



-120,0 



-160,0 



-240,0- 













Meson 

TT 


a 

P 

CJ 

■ ■ ■ ■ Tot 



0,0 1,0 

r [fm] 

(17): AcN(^Si) o E*iV("Di) 



2.0 



240,0 
200.0 
160.0 
^ 120.0 



0.0 
-40.0 





Meson 

7T 


a 

P 

— - - CJ 
■ ■ ■ ■ Tot 









0,0 1,0 

r [fm] 

(24): EciV(^Si) O AcNi^Di) 



2.0 



> 
CD 



210,0 
180,0 
150,0-1 
120,0 
90,0 
60,0 
30,0 
0,0 



p 

a 

Tot 



ao 1,0 

r [fm] 

(13): AcN^Si) ^ E*iV(=*5'i) 



2.0 



> 



15,0 
0,0 

-15,0 
-30,0 
-45,0 
-60,0 



/ "~ ^ ^ 










Meson 




7T 


CT 

P 

CJ 

Tot 



0,0 1,0 

r [fm] 

(16): AcN^Si) ^ E*iV(3£)i) 



2.0 



> 
0) 



1:^ 




1,0 

r [fm] 

(23): EciV(^Si) o T,*N{^Si) 



2,0 



> 




1,0 

r [fm] 

(25): E,7V(='Si) ^ T.cN{^Di) 



2.0 



10.0- 
0.0 
-10.0 
-20.0 
-30.0 
-40.0 
-50.0 















[MeV] 










Meson 






7T 






CT 

P 

CJ 

■ ■ ■ ■ Tot 





0.0 1.0 

r [fm] 

(26): E.iVf Si) o E*iV(=*Di) 



2,0 



15.0 
0.0 
-15.0 
-30.0 
-45.0 
-60.0 













Meson 




7T 


a 

P 

u 

Tot 



0.0 1.0 

r [fm] 

(34): E*iV(='Si) ^ AcN^Di) 



2.0 



75.0 
60.0 
45.0 
30.0 
15.0 
0.0 
-15.0 
-30.0 





Meson 

7T 


a 

P 

— - - CJ 

■ ■ ■ ■ Tot 






\ ^ ■' 

\ 



0.0 1.0 

r [fm] 

(36): E*iV(35'i) ^ E*Ar(3Di) 



2.0 




1.0 

r [fm] 

(45): AciV(=*Di) o EcAf(^£>i) 



2.0 



0.0 
-30.0 
-60.0 
-90.0 
-120.0 
-150.0 



/ ^ 










Meson 

7T 


a 

P 

CJ 

■ ■ ■ ■ Tot 



0.0 1.0 

r [fm] 

(27): EciV(^S'i) ^ T,*NCDi) 



2.0 



> 
CD 




1.0 

r [fm] 

(35): E*iV(='5'i) ^ EciV(^L»i) 



> 
0) 




1.0 

r [fm] 

(37): T.*N{^Si) o E*iV(5Di) 



2.0 



> 



315.0- 
270.0 
225.0 
180.0 
135.0 
90.0 
45.0 



0.0 



p 

CJ 

Tot 



0.0 1.0 

r [fm] 

(46): AcNi^Di) o E*Af(='Di) 



2.0 




2,0 



(47): AcN('^Di) ^ ElN^Di) 




2.0 



> 
0) 




0.0 1.0 

r [fm] 

(56): E,iV(^Di) o T.*NCDi) 



2.0 



> 




(57): T.cN(:'Di) ^ T.*N{-'Di 



1.0 

r [fm] 

(67): E*iV(^Di) « E*iV(5Di) 



FIG. 5: The potentials of different channels for the J =1^ case with — Aa 



Av 



1 GeV. (cont.) 



24 





GeV 


(B.E. in MeV) 


JP= 1 + 


GeV 


(B.E. 


in MeV) 




















^Do 






3Sl,3Dl,5Dl 






0.47 

10.27 




ScN 3.4 






ScN 3.4 








33.66 
(w/o) 


'So 






3Si, 3Di 


0.08 ■ 


4.09 

17.34' 


44.69 '■, 




3.3 






3.3 






91.52 


(w/o) 


A N 
















'So 

3.2 


0.64 ■ 


6.16 ■-. 

18.51 

38.88 

68.29 \ 


^Si, (3Di) 

3.2 




1.17 

20.64'' 


7239" 




3.1 




107.64 
(w/) 


3.1 











1.2 1.3 1.4 1.5 1.6 1.7 1.0 1.2 1.4 1.6 1.8 2.0 

(GeV) (GeV) 

= 0+ = 1+ 



FIG. 6: The sensitivity of tiie binding energy (B.E.) to the cutoff in the OPEP model for the = 0"'" and — 1+ cases. 
The cases without (w/o) and with (w/) channel coupling are both shown. (^Di) means there is no S — D mixing when one 
considers only the AcA^ channel. 



jp=o+ 


GeV 


(B.E. in MeV) 






JP= 1+ 


GeV 




(B.E. in MeV) 




























5Do 

EcN 3.4 










3Sl,3Dl,5Dl 
EcN 3.4 






2005-270^ 


42.89- 


,55.02 


'So 




2.22 •-. , 

14.22--, 

'3T35- 


52158- 




3Si,3Di 


2T84" 


2630'' 






7736 
(w/o) 


3.3 








74.62 
(w/o) 


3.3 




783r. 








AcN 


1.24 








AcN 




1.24 


''■151-59 (w/o) 






'So 

3.2 


0.12 

1330', 


2T()7-,, 

'5235, 


6845- 




^Si, (3Di) 

3.2 




4730'. 


4635-, 

6845>. 






3.1 




l23T4(*'> 




'5138 
(w/o) 


3.1 






, „ 9T38 

106716 


(w/o) 






0.8 0.9 


1.0 1.1 1.2 

Acom (GeV) 


1.3 


1.4 




0.8 


0.9 1.0 


1.1 1.2 1.3 1.4 
Acom (GeV) 


1.5 


1.6 1.7 



= 0+ = 1+ 

FIG. 7: The sensitivity of the binding energy (B.E.) to the cutoff A,r = Ao- = Avcc = Acom in the OBEP model for the = O"*" 
and = 1+ cases. The cases without (w/o) and with (w/) channel coupling are both shown. (^Di) means there is no S — D 
mixing when one considers only the AcA channel. 



25 



J''=0+ |GeV 



(B.E. in MeV) 



5Do 

EcN3.4- 



21.28 



'So 



017 "41)6- ■ 



ITTV 



44(1-. 



3.3- 



(w/o) 

5735-.':71.84 

145:97 
(w/o) 



AcN 



0.12 



3.2- 11.1^,^ 



8.80 • 



43lB> 



5T59-. 



79.43 
(w/o) 



3Sl,3Dl,5Dl 

EcN 3.4 



3Si, 3D 



3.3 



3Si, (3Di) 
3.2 



3.1 



1.2 1.5 1.8 2.1 2.4 2.7 3.0 3.3 3.6 3.9 

a 



GeV 



(B.E. in MeV) 



6.39 ■, 



'-.(w/o) 
38.38 



2.46 -, 



5T3<5- 



95777' 



0.65 150.55 



(w/o) 



i2:26-, 1T33-V 
4a9S- 



30M- 



(w/) 



69.86 
(w/o) 



1.0 1.3 1.6 1.9 2.2 2.5 2.8 3.1 3.4 3.7 

a 



FIG. 8: The sensitivity of the binding energy (B.E.) to the parameter a in the OBEP model for the = 0+ and = 1+ 
cases. The cases without (w/o) and with (w/) channel coupling are both shown. {^Di) means there is no S — D mixing when 
one considers only the AcN channel. 



